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La presente invention a pour objet une methode pour deformer 
graduellement des representations ou realisations, generes par simulation 
sequentielle, d'un modele stochastique non necessairement gaussien d'un 
milieu heterogene, en se basant sur un algorithme de deformation graduelle de 
modeles stochastiques gaussiens. 

La methode selon 1' invention trouve des applications dans la 
modelisation de zones souterraines oil il s'agit de generer des representations 
montrant comment est distribute une certaine grandeur physique dans une 
zone du sous-sol (la permeabilite z par exemple), qui soient compatibles au 
mieux avec des donnees observees ou mesurees : donnees geologiques, 
enregistrements sismiques, mesures obtenues dans des puits notamment 
mesures des variations au cours du temps de la pression et du debit de fluides 
issus d'un gisement etc. 

Etat de la technique 

Par la demande de brevet FR 98/09.018, on connait une methode pour 
deformer graduellement un modele stochastique (de type gaussien ou 
apparente) d'un milieu heterogene tel qu'une zone souterraine, contraint par 
une ensemble de parametres relatifs a la structure du milieu. Cette methode 
comporte le tirage d'un nombre p (p=2 par exemple) de realisations (ou 
representations) independantes du modele ou d'une partie au moins du modele 
choisi du milieu, parmi F ensemble de toutes les realisations possibles et une 
ou plusieurs etapes iteratives de deformation graduelle du modele en 
effectuant une ou plusieurs combinaisons lineaires successives de p 
realisations initiales independantes entre elles puis des realisations composites 
successivement obtenues eventuellement avec de nouveaux tirages, etc., les 




coefficients de cette combinaison etant tels que la somme de leurs carres est 
egale a 1. 

Les modeles gaussiens ou apparentes conviennent pour modeliser des 
champs de grandeurs continues et, de ce fait, sont mal adaptes a la 
5 modelisation de zones traversees par des reseaux de fractures ou des systemes 
de chenaux, par exemple. 

On connait par ailleurs les algorithmes dits de simulation sequehtielle, 
les plus largement utilises parmi les algorithmes de simulation geostatistique. 
Bien qu'ils conviennent particulierement pour la simulation des modeles 
10 gaussiens, il n'impliquent pas en principe une limitation a ce type de modele. 

On forme une representation geostatistique d'une zone soutetraine par 
exemple en la subdivisant par un reseau k N mailles et en determinant un 
vecteur aleatoire a N dimensions Z = (Zi, Z 2 ,.... Zk) correspondant au mieux a 
des mesures ou observations obtenues sur la zone. Comme le montre par 

15 exemple Johnson, M.E. ; in « Multivariate Statistical Simulation » ; Wiley & 
Sons, New York, 1987, cette approche reduit le probleme de la creation d ? un 
vecteur a N dimensions a une serie de N problemes mono-dimensionnels. Un 
tel vecteur aleatoire n ? est ni necessairement multi-gaussien ni stationnaire. La 
simulation sequentielle de Z implique d'abord la definition d'un ordre selon 

20 lequel les N elements (Zi, Zo,.... Zn) du vecteur Z sont generes Tun apres 
T autre. Sans rien perdre de generalite, on suppose que les N elements de Z 
sont generes en sequence depuis Zi jusqu'a Zn. Pour tirer une valeur de chaque 
element Zi, (i = 1, .... N), il faut realiser les operations suivantes : 



a) construire la distribution de Zi conditionnee par (Z !? Z2.. Z[. { ) 



F c (z i ) = P(Z i <z i /Z ly Z 2 



Z M ) ' et 



b) tirer une valeur de Z L de la distribution F c (zl). 

Dans la pratique de la geostatistique, la simulation sequentielle est 
utilisee firequemment pour engendrer des vecteurs multi-gaussiens et des 
vecteurs d'indicateur non-gaussiens. La tache principale de la simulation 
sequentielle est la determination des distributions conditionnelles F c (zO (i = 1, 
N). Des algorithmes et logiciels pour estimer ces distributions sont decrits 
par exemple dans : 

- Deutsch, C.V. et al, « GSLEB (Geostatistical Software Library) and User's 
Guide » ; Oxford University Press, New York, Oxford 1992. 

En ce qui conceme le tirage des valeurs a partir de la distribution F c (zO, 
il existe aussi un vaste ensemble d' algorithmes connus. 

On considere la methode de distribution inverse par laquelle on obtient 
une realisation de Z^: z i = ou ui est tire d'une distribution uniforme 

entre 0 et 1. II s'ensuit qu'une realisation du vecteur Z correspond a une 
realisation du vecteur U dont les elements Ui ? U2, ....,Un. sont mutuellement 
independants et uniformement distribues entre 0 et 1 . 

On s'aperqoit qu'une simulation sequentielle est une operation S qui 
transforme un vecteur uniforme U= (U l? U2, ....,Un) en un vecteur structure Z 



= (Z 1? Z 2 



Z N ): 



Z = S(U) 



(1) 



Le probleme que represente la contrainte d'un vecteur Z a differents 
types de donnees peut etre resolu en contraignant les distributions 



• 
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conditionnelles F c (zO (i = 1, N) et/ou le vecteur uniforme U= (Ui, U2, 



Des travaux recents sur ralgorithme sequentiel ont ete concentres sur 
T amelioration de 1' estimation des distributions conditionnelles F c (zO par 
5 integration de donnees geologiques et de donnees sismiques. 

On peut citer par exemple le papier de Zhu, H. et al ; « Formatting and 
Integrating Soft Data : Stochastic Imaging via the Markov-Bayes Algorithm » 
in Soares, A., Ed. Geostatistics Troia 92, vol.1 : Kluwer Acad. Publ., 
Dordrecht, The Netherlands, p.1-12, 1993. 

10 Cependant, cette approche ne peut etre etendue a F integration de 

donnees non lineaires comme les pressions issues des tests de puits et les 
historiques de production, sauf h imposer une severe linearisation. En outre, 
comme toute combinaison de vecteurs uniformes U ne donne pas un vecteur 
uniforme, la methode de deformation graduelle d'un modele stochastique 

15 developpee dans le brevet deja cite, ne peut etre directement appliquee dans le 
cadre de la technique sequentielle rappelee ci-dessus. 

Aussi, la methode selon T invention permet de rendre compatible les 
deux approches i.e. d'etendre le formalisme developpe dans le brevet 
precedemment cite, a la deformation graduelle de realisations, generees par 
20 simulation sequentielle, d'un modele stochastique non necessairement 
gaussien. 



■U N ). 



Definition de la methode 



La methode permet de deformer graduellement une representation ou 
realisation, generee par simulation sequentielle. d un modele stochastique non 




necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre a un ensemble de donnees 
collectees dans le milieu par des mesures et observations prealables, relatives a 
son etat ou sa structure. 

5 Elle est caracterisee en ce qu'elle comporte 1' application d'un 

algorithme de deformation gradueile d'un modele stochastique a un vecteur 
gaussien (Y) ayant un nombre N de variables mutueilement independantes qui 
est relie a un vecteur uniforme (U) a N variables uniformes mutueilement 
independantes par la fonction de repartition gaussienne (G), de fagon k definir 

io une chaine de realisations u(t) du vecteur (U), et T utilisation de ces 
realisations u(t) pour generer des realisations z(t) de cette grandeur physique 
que Ton cale par rapport aux donnees (non lineaires). 

Suivant un premier mode de mise en ceuvre, on definit la chaine de 
realisations u(t) du vecteur uniforme (U) a partir d'une combinaison lineaire 
15 de realisations du vecteur gaussien (Y) dont les coefficients de combinaison 
sont tels que la somme de leurs carres est egale a 1' unite. 

Suivant un autre mode de mise en oeuvre, on realise une deformation 
gradueile d ? un nombre n de parties du modele representatif du milieu 
heterogene tout en preservant la continuite entre ces n parties du modele, par 
20 une subdivision du vecteur uniforme (U) en un nombre n de sous-vecteurs 
mutueilement independants. 

D'autres caracteristiques et avantages de la methode selon rinvention, 
apparaitront a la lecture de la description ci-apres d*un exemple non limitatif 
de realisation, en se referant aux dessins annexes ou : 



- la Fig.l montre la couche mediane d'une realisation d'un modele de fades 
servant de reference, generee par simulation sequentielle d'indicatrices ; 

- la Fig. 2 montre la variation en fonction du temps de la pression obtenue 
dans un test de puits du modele de reference ; 

5 - les Fig.3A a 3E montrent respectivement cinq realisations initiales de la 
couche mediane d'une zone reservoir, contraintes seulement par les facies le 
long du puit ; 

- les Fig. 4A a 4E montrent respectivement pour ces cinq realisations, les 
courbes des pressions au puits dans le modele de reference comparees h 

10 celles correspondant aux modeles initiaux ; 

- les Fig.5A k 5E montrent respectivement cinq realisations de la couche 
mediane du modele de facies conditionnees aux facies le long du puits et 
calees par rapport a la courbe de pression obtenue par test de puits ; 

~ les Fig. 6A a 6E montrent respectivement pour les cinq realisations, les 
15 courbes des pressions au puits dans le modele de reference comparees a 
celles correspondant aux modeles cales ; 

- les Fig.7A a 7E montrent respectivement comment les fonctions objectif 
correspondant respectivement a ces cinq exemples, varient avec le nombre 
d' iterations ; 

20 - les Fig. 8A a 8E montrent les deformations graduelles operees par un 
changement de coefficient d' anisotropic sur un modele a trois facies genere 
par simulation sequentielle d'indicatrices; et 
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- les Fig. 9 A a 9E montrent les deformations graduelles locales d'un modele a 
trois facies, generees par simulation sequentielle d' indie atrices. 

Description detaillee de la methode 

On considere une zone d' etude que Ton subdivise par tine grille a N 
mailles et Ton cherche a construire des realisations ou representations d'un 
modele stochastique d'une certaine grandeur physique z representant par 
exemple la permeabilite des formations dans la zone. Le modele recherche doit 
s'ajuster a des donnees obtenues par des mesures ou observations en un certain 
nombre de points, et notamment a des donnees non lineaires. 

Le calage d'un modele stochastique a des donnees non lineaires peut 
etre formulee comme un probleme d' optimisation. On designe par f° bs = (f t obs , 
f2° bs , f3° bs f P ° bs ), le vecteur des donnees non lineaires observees ou mesurees 

dans le milieu etudie (la zone reservoir), et par f = (f! , f 2 , f 3 , f p ), le 

vecteur correspondant des reponses du modele stochastique de la permeabilite 
Z = (Zi, Z 2 ,..., Zn). Le probleme de contraindre le modele stochastique Z par 
les observations consiste a generer une realisation z de Z qui reduit a une 
valeur assez faible une fonction objectif que Ton definit comme la somme des 
erreurs quadratiques ponderees des reponses du modele par rapport aux 
observations ou mesures dans la zone-reservoir ? soit : 

ou toy represente le poids attribue a la reponse fj. Les fonctions fi (i=l. 2 ? 

p) ainsi que la fonction objectif O sont des fonctions du vecteur Z. On se 
trouve devant un probleme d' optimisation de dimension N. 
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Pour etendre le formalisme developpe dans le brevet deja cite, a la 
deformation graduelle de realisations generees par simulation sequentielle non 
necessairement gaussienne, on va partir d'un vecteur gaussien a N variables Yi 
avec i = 1, 2, N, mutuellement independantes, de moyenne nulle et de 
5 variance egale a 1' unite, et on va definir N variables uniformes mutuellement 
independantes Ui, U 2) U 3 ,...U N par : 

Ui = G(Y0 V i= 1,2, N 

ou G designe la fonction de repartition gaussienne standardisee. 

Ceci pose, on va appliquer ralgorithme de deformation graduelle 

io developpe dans un cadre gaussien, au vecteur gaussien Y = (Yi, Y 2 , Y N ), 

pour construire une chaine continue de realisations du vecteur uniforme U = 

(Ui, U 2 U N ). Etant donne deux realisations independantes y a et y b de Y, on 

definit la chaine de realisations u(t) du vecteur U que l'on obtient selon la 
relation : 

15 u(0 = G( y a cosf + y b sin t ) (2) 

Pour chaque t, u(t) est une realisation du vecteur U. Ensuite, par 

echantillonnage de la distribution conditionnelle F f (s.)(i=l, 2, N) en 

utilisant les elements du vecteur u(t), on obtient un vecteur z(t) qui est. pour 
chaque t, une realisation du vecteur aleatoire Z. En consequence, on peut 
20 ajuster le parametre t comme dans le cas gaussien de fa?on a caler z(t) a des 
donnees non lineaires. Cette procedure est iteree jusqu'a atteindre un calage 
satisfaisant. 




Calage d'un modele de facies a des donnees de pression obtenues par 
des tests de puits 

Pour illustrer 1' application de la methode d ? optimisation stochastique 
ci-dessus definie, on va chercher a caler un modele stochastique de reservoir k 
des donnees de pressions obtenues par des tests de puits. La construction du 
modele de reservoir est derivee d'une formation petroliere reelle qui comporte 
trois types de facies : deux facies reservoir de bonne qualite (facies 1 et 2) et 
un type de facies reservoir de tres mauvaise qualite (facies 3). La table 1 
definit les proprietes petrophysiques des trois facies : 





K x (md) 


K y (md) 


K.(md) 


<£(%) 


c,(\Q- s bar- 1 ) 


facies 1 


10. 


10. 


10. 


17 


2.1857 


facies 2 


1. 


1. 


1. 


14 


2.0003 


facies 3 


0.1 


0.1 


0.001 


9 


1.8148 



Pour representer la distribution specifique de facies de la formation 
petrolifere, on genere d'abord une realisation binaire pour representer le facies 
3 et son complement. Ensuite, dans la partie complementaire du facies 3, on 
genere une autre realisation binaire independante de la premiere pour 
representer les facies 1 et 2. La formation est discretisee par un maillage 
regulier de 60x59x15 blocs de taille 15mxl5mxl.5m. Un modele de 
variogramme exponentiel est utilise pour estimer les distributions 
conditionnelles. La direction d'anisotropie principale est diagonale par rapport 
au maillage. Les portees du variogramme du facies 3 selon les trois directions 




d'anisotropie sont respectivement de 300m, 80m et 3m. Les portees du 
variogramme des facies 1 et 2 selon les trois directions d'anisotropie sont 
respectivement de 150m, 40m et 1.5m. Les proportions des facies 1, 2, 3 sont 
respectivement de 6%, 16% et 78%. 

5 On a realise un test de puits au moyen d'un simulateur de test de puits k 

differences finies, de la maniere decrite par : 

Blanc, G. et al : « Building Geostatistical Models Constrained by Dynamic 
Data - A Posteriori Constraints » in SPE 35478, Proc. NPF/SPE European 3D 
Reservoir Modelling Conference, Stavenger, Norway, 1996. 

10 On voit sur la Fig.l la couche mediane d'une realisation qui sert de 

modele de reference pour cette validation. La section du puits qui a ete 
perforee traverse horizontalement la couche mediane du modele de reservoir le 
long de l'axe x. Le diametre du puits est de 7.85cm, la capacite du puits est 
nulle et les facteurs parietaux (skin) des facies 1, 2 et 3 sont respectivement 0, 

15 3 et 50. Le test de puits synthetique dure 240 jours avec un debit constant de 5 
m3/jour de fagon k investiguer presque la totalite du champ petrolifere. La 
Fig. 2 montre la variation de pression en fonction du temps. 

Le but etait de construire des realisations du modele de facies contraint 
par les facies rencontre le long du puits et par la courbe de pression obtenue 

20 lors du test de puits. La fonction objectif est definie comme la somme des 
differences quadratiques entre les reponses en pression du modele de reference 
et les reponses en pression de la realisation. Du fait que le compottement 
dynamique du modele de reservoir est principalement controle par le contraste 
entre les facies reservoir de bonne et de mauvaise qualites, la realisation 

25 binaire utilisee pour generer les facies 1 et 2 a ete fixee d'abord et seule la 



realisation binaire utilisee pour generer le fades 3 a ete deformee pour le 
calage des donnees de pression. 

Les reponses en pression resultant des tests de puits pour les cinq 
realisations des Fig.3A a 3E, sont differentes de celle du modele de reference, 

5 comme le montrent les Fig. 4A a 4E. Partant respectivement de ces 5 
realisations independantes, en utilisant la methdde de calage iteratif ci-dessus, 
on obtient apres quelques iterations, cinq realisations calees (Fig. 5 A a 5E) 
pour lesqueiles les courbes de pression correspondantes sont tout a fait en 
accord avec celles du modele de reference, comme le montrent les Fig.6A & 

io 6E. 

Deformation graduelle par rapport aux parametres structuraux 

Dans de nombreux cas, on ne dispose pas de donnees suffisantes pour 
en deduire les parametres structuraux du modele stochastique : moyenne, 
variance, fonction de covariance, etc. Ces parametres structuraux sont souvent 

15 donnes en termes d'intervalles ou de distributions a priori. Si leurs valeurs sont 
fausses, il est vain de rechercher une realisation calee sur des donnees non 
lineaires. Aussi, il est essentiel pour les applications de pouvoir proceder a une 
deformation graduelle d'une realisation avec modification simultanee de 
nombres aleatoires et des parametres structuraux. L ? algorithme de simulation 

20 sequentielle definie par T equation (1) fournit cette possibility de changer 
simultanement ou separement Toperateur structurel S et le vecteur uniforme U. 
Les Fig.8A-8E montrent les deformations graduelles que Ton obtient pour une 
realisation fixe du vecteur uniforme U quand on change le coefficient 
d' anisotropic 
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Deformation graduelle locale ou regionalisee 

Quand les observations sont reparties sur des zones difKrentes d'une 
formation etudiee, un calage qui utiliserait une deformation globale risquerait 
d'etre inefficace car T amelioration de 1' accord que Ton obtiendrait sur une 
zone pourrait le deteriorer sur une autre zone. D est done bien preferable 
d'appliquer des deformations graduelles zone par zone. On considere une 
subdivision du vecteur U en un certain nombre n de sous-vecteurs U 1 , U 2 ? .... ? 
U n mutuellement independants ce qui rend possible de realiser 
individuellement leur deformation graduelle. En appliquant separement a 
chacun des sous-vecteurs U 1 , U 2 ,...., U n , Talgorithme de deformation 
graduelle, on obtient une fonction de dimension n, du vecteur uniforme U : 

G(Y* cos t 2 + Y* sin t 2 ) 
G(Y a n cost n + Y b "smt n )_ 

ou Y l a et Y£ pour tout i =1, 2, n, sont des sous-vecteurs gaussiens 
independants. Pour un ensemble de realisations de et Yf , on resout un 
probleme d' optimisation de n parametres ti, t 2 , t n pour obtenir une 
realisation qui ameliore, sinon maintient, le calage des donnees. Cette 
procedure peut etre iteree jusqu'k atteindre un calage satisfaisant. 

En realisant des deformations locales graduelles, on ameliore de fa?on 
significative la vitesse de calage dans tous les cas ou les mesures ou 
observations sont reparties dans differentes zones du milieu. 



U\t x ) 
U z (t z ) 

U n (t n ) 




On voit bien Feffet de cette deformation locale graduelle sur le modele 
a trois facies, des Fig. 9A a 9E ou seule la partie inferieure gauche delimitee 
est affectee. 

La methode seloa 1' invention peut etre facilement generalisee a la 
deformation graduelle d'une representation ou realisation d'un modele 
stochastique quelconque, car generer une realisation d'un tel mqdele 
stochastique se ramene toujours a la generation des nombres uniformes. 
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REVENDICATIONS 



• 1) Methode pour deformer graduellement une representation ou 
realisation, generee par simulation sequentielle, d'un modele stochastique non 
necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, arm de la contraindre a un ensemble de dbnnees 
collectees dans le milieu par des mesures et observations prdalables, relatives a 
son etat ou sa structure, caracterisee en ce qu'elle comporte 1' application d'un 
algorithme de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) a N variables mutuellement independantes qui est relie a un 
vecteur uniforme U a N variables uniformes mutuellement independantes par 
la fonction de repartition gaussienne (G), de facon a construire une chatne de 
realisations u(t) du vecteur U, et l'utilisation de ces realisations u(t) pour 
generer des realisations z(t) de cette grandeur physique, que l'on cale par 
rapport aux donnees. 

o 2) Methode selon la revendication 1, caracterisee en ce que Ton definit 
une chaine de realisations u(t) du vecteur (U) a partir d'une combinaison 
lineaire de realisations du vecteur gaussien (Y) dont les coefficients de 
combinaison sont tels que la somme de leurs carres est egale a 1' unite. 

* 3) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle du modele representatif du milieu 
heterogene, simultanement par rapport aux parametres structuraux et aux 
nombres aleatoires. 



# 



c 4) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle separee d'un nombre n de parties du 
modele representatif du milieu heterogene tout en preservant la continuite 
entre ces n parties du modele, par une subdivision du vecteur uniforme en n 
sous-vecteurs mutuellement independants. 




Fig.8A Fig.8B Fig.8C Fig.8D Fig.8E 




Fig.9A Fig.9B Fig.9C Fig.9D Fig.9E 



